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Abstract. The Color Glass Condensate is a theory of the dynamical properties of partons in the Regge limit 
of QCD: xsj — > 0, Q 2 >> Aq CD — fixed and the center of mass energy squared s — > oo. We provide a brief 
introduction to the theoretical ideas underlying the Color Glass Condensate and discuss the application of 
these ideas to high energy scattering in QCD. 
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1 Introduction 

The study of the properties of the strong interactions 
in the asymptotic Bjorken limit of momentum transfer 
squared Q 2 — > oo, the center of mass energy squared 
s — > oo, and the Bjorken variable xb} ~ Q 2 /s — fixed 
has proved to be one of the most creative ideas in theo- 
retical physics £Q. Relatively little work has been done in 
the other high energy limit, namely, XBj — > 0, s — > oo and 
Q 2 = fixed. This limit of the strong interactions, which 
we shall call the Regge limit, was studied intensively in 
the 60's and indeed led eventually to string theory. The 
reason these studies fell into disfavor in the strong inter- 
actions was that there was no small parameter in these 
studies (in modern parlance, Q 2 < ^qcd)- 

With the advent of the collider era, we can now probe 
a wide window of physics where s » Q 2 >> ^1q CD - 
In fact, this "window" describes the bulk of the high en- 
ergy cross-section. One therefore has finally the possibility 
of studying the properties of the Regge limit of the the- 
ory using weak coupling methods. In this limit, the hadron 
behaves like matter that's dense but weakly coupled-not 
dissimilar to much of condensed matter physics |2B2]- 

In Regge asymptotics, the number of partons increases 
rapidly due to QCD bremsstrahlung. This growth is de- 
scribed, in the leading logarithmic approximation in x, by 
the BFKL equation fj] . Since the typical size of the par- 
tons in this limit is of order 1/Q 2 , the hadron becomes 
closely packed when the number of partons is of order 
R 2 Q 2 . In fact, this corresponds to an occupation number 
/ ~ 1/ois- When the density of partons is of this order, re- 
pulsive many body "recombination" and screening effects 
compete with QCD Bremsstrahlung leading to a satura- 
tion of the number of partons in the hadron's wavefunc- 
tion |S]n>|I7]. The saturation of partons of different sizes 
happens at different values of x. The scale at which this 
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occurs is the saturation scale Q s (x)-a, dynamically gener- 
ated semi-hard scale that controls the dynamics of physics 
in this regime of QCD. 

In the language of the Operator Product Regime (OPE), 
the line Q = Q s (x) in the x-Q 2 plane denotes the regime 
beyond which (when approached from high Q 2 ) higher 
twist effects become important. Recall that the OPE is 
best formulated in the Bj-limit where higher twists are 
power suppressed and can be forgotten. The opposite is 
true in the Regge limit. Since the number of twist opera- 
tors grows (nearly) exponentially with the twist, the OPE 
quickly becomes unwieldy. Thus to describe physics in this 
regime we need a new organizing principle in QCD beyond 
the OPE. 



2 A classical effective theory for high energy 
QCD 

A way out was suggested when it was realized that the 
physics of high parton densities could be formulated as a 
classical effective theory When a quantum field the- 
ory is formulated on the light cone, one realizes that there 
is a formal Born-Oppenheimer separation between large 
x and small x modes 9 which are respectively the slow 
and fast modes in the effective theory. Thus on the time 
scale of the "wee" parton small x fields, the large x par- 
tons can be viewed as static charges. Since these are color 
charges, they cannot be integrated out of the theory but 
must be viewed as sources of color charge for the dynami- 
cal wee fields. With this dynamical principle in mind, one 
can write down an effective action for wee partons in QCD 
at high energies. The generating functional of wee partons 
has the form 



Z[j}= [dp]W A+ [p] 



F [dA]S(A+)e 



iS[A,p]-j-A 



J A+ [dA]S(A+)e iS ^ 



(1) 
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where the wee parton action has the form 



S[A,p} = 



-1 



d 2 x±dx S(x ) 



x Ti(p(x ± )U-oo,oo{A-}) 



(2) 



In Eq. ^ P is a classical color charge density (more on 
this shortly) of the static sources and W[p] is a weight 
functional of sources (which sit at momenta k + > A + : 
note, x = fc + /-Phadron)- The sources are coupled to the 
dynamical wee gluon fields (which in turn sit at k + < 
A + ) via the gauge invariant term 1 which is the first term 
on the RHS of Eq. [3 The second term in Eq. H is the 
QCD field strength tensor squared-thus the wee gluons 
are treated in full generality in this effective theory, which 
is formulated in the light cone gauge A + = 0. The source 
j is an external source-derivatives taken with respect to 
this source (with the source then put to zero) generate 
correlation functions in the usual fashion. 

We have not justified thus far why the sources are clas- 
sical. The argument for this is subtle and follows from 
a coarse graining of the effective action to only include 
modes of interest. For large nuclei, or at small x, the wee 
partons couple to a large number of sources. For a large 
nucleus, it can be shown explicitly that this source den- 
sity is classical Further, it was conjectured that the 
weight functional for a large nucleus was a Gaussian in the 
source density (corresponding to the quadratic Casimir 
operator) [HUE]- This was shown explicitly recently to be 
the correct- albeit with small corrections for SU(-/V C ) com- 
ing from the N c — 2 higher Casimir operators 

For a large nucleus, the variance of the Gaussian, the 
color charge squared per unit area p\, proportional to 
A 1 / 3 , is a large scale-and is the only scale in the effective 
action 2 . Thus for p 2 A » ^1q CD , «s(^) << 1, and one 
can compute the properties of the theory in Eq. [2 in weak 
coupling. 

By evaluating the saddle point of the action in Eq. [2J 
one can compute the classical distribution of gluons in the 
nucleus. The Yang-Mills equations can be solved analyti- 
cally to obtain the classical field of the nucleus as a func- 
tion of p: A c i (p) f5fr2yr3j . From the generating functional 
in Eq. ^ one obtains for the two point correlator, 



< AA >= / [dp] W A + [p] A cL (p) A d .{p) . 



(3) 



function and z = k\/Q 2 .) A very nice expression for the 
classical field of the nucleus containing these two lim- 
its was presented by Triantafyllopoulos at this confer- 
ence [T3] . 

We are now in a position to discuss why a high en- 
ergy hadron behaves like a Color Glass Condensate 0- 
The "color" is obvious since the degrees of freedom, the 
partons, are colored. It is a glass because the stochastic 
sources (frozen on time scales much larger than the wee 
parton time scales) induce a stochastic (space-time depen- 
dent) coupling between the partons under quantum evolu- 
tion (to be discussed in the next section)-this is analogous 
to a spin glass. Finally, the matter is a condensate since 
the wee partons have large occupation numbers (of order 
1/as) and have momenta peaked about Q s . As we will 
discuss, these properties are enhanced by quantum evolu- 
tion in x. The classical field retains its structure-while the 
saturation scale grows: Q s (x') > Q s (x) for x' < x. 



3 Quantum evolution a la JIMWLK and BK 

Small fluctuations about the effective action in Eq.|2|were 
first considered in Ref. [T5|. It was discovered that these 
gave large corrections of order ln(I/:r). In particular, 
this suggested that the Gaussian weight functional was 
fragile under quantum evolution of the sources 3 . A Wilso- 
nian renormalization group (RG) approach was developed 
to systematically treat these corrections ^j. The basic in- 
gredients of this approach are as follows. Begin with the 
generating functional in Eq. ^ &t some A + , with an ini- 
tial source distribution W[/>]. Perform small fluctuations 
about the classical saddle point of the effective action, in- 
tegrating out momentum modes in the region A l+ < k + < 
A + , ensuring that A' + is such that as \n(A + /A' + ) << I. 
The action reproduces itself at the new scale A , al- 
beit with a charge density p' = p + Sp, and Wa+ [p] — ► 
W^/+ [p']. The change of the weight functional W[p] with x 
is described by the JIMWLK- non-linear RG equation \IQ 
which we shall not write explicitly here. 

The JIMWLK equations form an infinite hierarchy (anal- 
ogous to the BBGKY hierarchy in statistical mechanics) 
of ordinary differential equations for the gluon correlators 
< A\A2 - ■ ■ A n >y, where Y = ln(I/x) is the rapidity. 
The expectation value of an operator O is defined to be 



From this expression one can determine, for Gaussian sources, 
the occupation number <ft — dN/nR 2 jdk\dy of wee par- 
tons in the classical field of the nucleus. One finds that 



[da}0[a}W Y [a] 



(4) 



for k± >> Ql, one has the Weizsacker- Williams spectrum 
4> ~ Q 2 s /k\ : while for k± < Q s , one has a complete resum- 
mation to all orders in k±, which gives <j> ~ -^ln(Q s /k±). 
(The behavior at low fcj_ can, more accurately, be repre- 
sented as ^--T(0, z) where r is the incomplete Gamma 

1 This is not the only possible gauge invariant coupling. An 
alternative form is given in Ref. [10 -it can be shown to repro- 
duce BFKL more efficiently. 

2 p\ is simply related in the classical theory to the saturation 
scale Qg via the relation Q 2 a = asN c p?A WQs/^qcd) 



where a — -^r p. The corresponding JIMWLK equation 
for this operator is 



d < Q[a] > Y 



1 

< 2 



Sa Y (y A 



<yx Sa Y (x±) 
-OH > Y ■ 



ab 

Ax ± ,y ± 



(5) 



We will return to this point in our discussion of the Cronin 
effect in Deuteron-Gold collisions at RHIC. 
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X here is a non-local object expressed in terms of path 
ordered (in rapidity) Wilson lines of a • This equation 
is analogous to a (generalized) functional Fokker- Planck 
equation, where Y is the " time" and \ is a generalized dif- 
fusion coefficient. This equation illustrates the stochastic 
properties of operators in the space of gauge fields at high 
energies. For the gluon density, which is proportional to 
a two-point function < a a (x±)a b (y±) >, one recovers the 
BFKL equation in the limit of low parton densities. 

As mentioned, the JIMWLK equations are master equa- 
tions for n-point correlators. Two point correlators of Wil- 
son lines are proportional to 4-point correlators and so 
on. The theory is conformal so it is not inconceivable 
that it is exactly solvable but this has not been done 
thus far. Preliminary numerical solutions have been ob- 
tained recently |17| but much work remains in that direc- 
tion. There is a mean field solution deep in the saturation 
regime |18| where one can show that the weight functional 
is a Gaussian-albeit a non local Gaussian with a variance 
proportional to kj_ for k\ < Q 2 . 

In the limit of large N c and large A (a^A 1 ^ >> 1), 
one can show that the hierarchy closes for the two point 
correlator of Wilson lines since the expectation value of 
the product of traces of Wilson lines factorizes into the 
product of the expectation values of the traces: 

< TT(V x V^ z )Tr{V z V^) >^< Tr(V x V}) > < Tr(V z Vj) > , 

(6) 

where x — "Pexp(J dz a a (z ,a^)T a ). Here V denotes 
path ordering in x~ and T a is the SU(3) generator in 
the adjoint representation. In Mueller's dipole picture 4 , 
the cross-section for a dipole scattering off a target can be 
expressed in terms of these 2-point dipole operators as [TH| 

a q q N {x,r±) = 2 J d 2 b Af Y {x,r±,b) , (7) 

where Ay, the imaginary part of the forward scattering 
amplitude, is defined to be Afy = 1 - jj- < Tr(V x Vj) >y. 
Note that the size of the dipole, r± = x± — y ± and 
b = (x± + y±)/2. The JIMWLK equation for the two 
point Wilson correlator is identical in the large A, large 
N c mean field limit to an equation derived independently 
by Balitsky and Kovchegov-the Balitsky-Kovchegov equa- 
tion |21j . which has the operator form 

^=a S /Cbfkl <8> {Ay - Afy } ■ (8) 

Here /Cbfkl is the well known BFKL kernel. When Af << 
1, the quadratic term is negligble and one has BFKL 
growth of the number of dipoles; when Af is close to unity, 
the growth saturates. The approach to unity can be com- 
puted analytically The B-K equation is the simplest 
equation including both the Bremsstrahlung responsible 
for the rapid growth of amplitudes at small x as well as 
the repulsive many body effects that lead to a saturation 
of this growth. 

4 See also Ref. EOT 



A saturation condition which fixes the amplitude at 
which this change in behavior is significant, say Af = 
1/2, determines the saturation scale. One obtains Q 2 S = 
QQexp(AF), where A = cas with c w 4.8. The saturation 
condition affects the overall normalization of this scale but 
does not affect the power A. In fixed coupling, the power 
A is large and there are large pre-asymptotic corrections 
to this relation-which die off only slowly as a function of 
Y. BFKL running coupling effects change the behavior of 
the saturation scale completely-one goes smoothly at large 
Y to Q 2 = Qq cxp(y / 2b c(Y + Y )) where b is the coef- 
ficient of the one-loop QCD /3-function. The state of the 
art computation of Q s is the work of Triantafyllopoulos, 
who obtained Q s by solving NLO-resummed BFKL in the 
presence of an absorptive boundary (which corresponds 
to the CGC) The pre-asymptotic effects are much 

smaller in this case and the coefficient A « 0.25 is very 
close to the value extracted from saturation model fits to 
the HERA data j22- 

No analytical solution of the BK equation in the en- 
tire kinematic region but there have been several numeri- 
cal studies at both fixed and running coupling |2*oTl25ll27j . 
These studies suggest that the solutions have a soliton like 
structure and that the saturation scale has the behavior 
discussed here. Geometrical scaling of solutions is seen for 
a wide window in rapidities. Running coupling effects, as 
suggested, are important and make the results of the com- 
putations more physically plausible. 

The soliton like structure is no accident, as was dis- 
covered by Munier and Peschanski who noticed that 
the BK-equation, in a diffusion approximation, bore a for- 
mal analogy to the FKPP equation describing the prop- 
agation of unstable non-linear wavefronts |3U| . In addi- 
tion, the full BK-equation lies in the universality class of 
the FKPP equation enabling one to extract the univer- 
sal properties of these equations (for instance the leading 
pre-asymptotic terms in the expression for the saturation 
scale). The power of this analogy was made manifest j2H] 
when it was realized that a stochastic generalization of 
the FKPP cquation-the sFKPP equation-could provide in- 
sights into impact parameter dependent fluctuations |31| 
in high energy QCD beyond the BK-equation. This is 
a very active area of research now, with several groups 
hunting for the Pomeron loops responsible for these fluc- 
tuations. The (rapidly evolving) state of the art of this 
subject is discussed in the talk by Edmond Iancu in these 
proceedings |52*]. 

To summarize, the Color Glass Condensate is a weak 
coupling effective theory describing the properties of hadron 
wavefunctions in QCD at high energies. Renormalization 
group equations-the JIMWLK equations-describe the be- 
havior of multi-parton correlations in the hadron wave- 
function as a function of rapidity. The theory has stochas- 
tic features closely analogous to the propagation of unsta- 
ble non-linear wave fronts in statistical mechanics. Recent 
work |32j is focused on trying to understand possible cor- 
rections beyond JIMWLK at low parton densities-which 
may be responsible for Pomeron loops. We now turn to 
the applications of this theory to hadronic scattering. 



4 



Raju Venugopalan: The Color Glass Condensate: an overview 



4 Hadronic scattering and k± 
the Color Glass Condensate 



factorization in 



Collinear factorization is the perturbative QCD mecha- 
nism to compute hard scattering. For instance, the cross- 
section in pp-collisions for di-jets with invariant mass M 2 ~ 
s >> Aq CI} , is a convolution of structure functions from 

each of the nucleons (evaluated at the scale M 2 ) times the 
probability that collinear partons with fcj^ = from the 
nuclei scatter to produce the di-jet. The structure func- 
tions are universal since they can be extracted from one set 
of experiments and input into another. Factorization theo- 
rems prove that this universality holds modulo power cor- 
rections in M 2 . At collider energies, a new window opens 
up where Aq CD << M 2 « s. In principle, cross-sections 
in this window can be computed in the collinear factor- 
ization language-however, one needs to sum up large log- 
arithmic corrections in s/M 2 . An alternative formalism is 
that of fcj_-factorization (331134) . where one has a convolu- 
tion of k± dependent "un-integrated" gluon distributions 
from the two hadrons with the hard scattering matrix. In 
this case, the in-coming partons from the wavefunctions 
have non-zero k±. It was suggested by Levin et al. 
that at high energies the typical k± is the saturation scale 
Q s . The rapidity dependence of the unintegrated distri- 
butions is given by the BFKL equation. However, unlike 
the structure functions, it has not been proven that these 
unintegrated distributions are universal functions. 

At small x, both the collinear factorization and k± fac- 
torization limits can be understood in a systematic way 
in the framework of the Color Glass Condensate. Rather 
than a convolution of probabilities, one has instead a col- 
lision of classical gauge fields. The expectation value of an 
operator O can be computed as 



< O >Y- 



[dpi] [d P2 ] W X1 [ Pl ] W X2 [p 2 ] 0( Pl , pa) , (9) 



where Y = 1ii(1/xf) and xp = x\ — x 2 . All operators at 
small x can be computed in the background classical field 
of the nucleus at small x. Quantum information, to lead- 
ing logarithms in x, is contained in the source functionals 
W X1 ( X2 ) [pi{p 2 )]- The operator O can be expressed in terms 
of gauge fields A^[px, p 2 ]{x). 

Inclusive gluon production in the CGC is computed 
by solving the Yang-Mills equations [D M , F» u ] a = J"' a , 
where 



J v = Pl 5{x-)5 U+ + p 2 5{x + )5 v 



(10) 



with initial conditions given by the Yang-Mills fields of 
the two nuclei before the collision. These are obtained 
self-consistently by matching the solutions of the Yang- 
Mills equations on the light cone |___]. The initial condi- 
tions are determined by requiring that singular terms in 
the matching vanish. Since we have argued in Section [3. 
that we can compute the Yang-Mills fields in the nuclei 
before the collision, the classical problem is in principle 
completely solvable. Quantum corrections not enhanced 
by powers of asln(l/x) can be included systematically. 



The terms so enhanced are absorbed into the weight func- 
tionals W[pi,2]- 

Hadronic scattering in the CGC can therefore be stud- 
ied through a systematic power counting in the density of 
sources in powers of p\. 2 lk\. x 2 - This power counting in 
fact is more relevant at high energies than whether the 
incoming projectile is a hadron or a nucleus. In addition, 
one can begin to study the applicability of both collinear 
and fcj_ factorization at small x in this approach. 



4.1 Gluon and quark production in the dilute/pp 
regime: (p pl / k\ p p2 / k\ « 1) 

The power counting here is applicable either to a proton at 
small x, or to a nucleus (whose parton density at high en- 
ergies is enhanced by A 1 / 3 ) at large transverse momenta. 
The relevant quantity here is Q Sl which, as one may re- 
call, is enhanced both for large A and small x. So as long 
as >> Q s » Aqcb, one can consider the proton or 
nucleus as being dilute. 

To lowest order in p v \jk\ and p P 2/k 2 _, one can com- 
pute inclusive gluon production analytically. This was first 
done in the A T = gauge |3fi| and subsequently in the 
Lorentz gauge d^A^ = 37 . At large transverse mo- 
menta, Q s << k±, the scattering can be expressed in 
a fcj_-factorized form. The inclusive cross-section is ex- 
pressed as the product of two unintegrated (k± dependent) 
distributions times the matrix element for the scatter- 
ing. The comparison of this result to the collinear pQCD 
gg — > gg process and the k± factorized gg — > g was per- 
formed in Ref. At this order, the result is equivalent 
to the perturbative QCD result first derived by Gunion 
and Bertsch jSHJ. This result for gluon production is sub- 
stantially modified, as we shall discuss shortly, by high 
parton density effects in the nuclei. 

k±_ factorization is a good assumption at large mo- 
menta for quark pair-production. This was worked out in 
the CGC approach by Francois Gelis and myself The 
result for inclusive quark pair production can be expressed 
in k i factorized form as 



da i 



d 2 ki± d 2 k 2 j 



dy p dy q d 2 p ± d 2 qi _ ~ J (2n) 2 (2^)^ (fclJ 

Tr (| m ab + ( fc l> fc 2 ;<?,£) | 



k 2 ± —p± — qj 



xipi(ki±)tp 2 (k 2 ±)- 



k 2 k 2 



(11) 



where 4>\ and <j) 2 are the unintegrated gluon distributions 
in the projectile and target respectively (with the gluon 

distribution defined as xG{x,Q 2 ) = L d(kj_) <p(x,k±)). 

The matrix element Tr (\ m -J ( kl , k 2 ; g,p)\ 2 ) is identical 

to the result derived in the k± -factorization approach |33l 

_____ , , Tr (|m" + (fci,fe 2 ;g.p)| 2 ) 

E3J. In the limit fen , k 2 ± -> 0, Vl ab J A 1 > is 
well defined-after integration over the azimuthal angles in 
Eq. rj_Q one obtains the usual matrix element \M\ 2 gg ^ q q : 
recovering the lowest order pQCD collinear factorization 
result. 
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4.2 Gluon and quark production in the 
semi-dense/pA region {p p /k\ « 1 p^/kj_ ~ 1). 

The power counting here is best applicable to asymmet- 
ric systems such as proton-nucleus collisions, which nat- 
urally satisfies the power counting for a wide range of 
energies. Of course, as one goes to extremely high en- 
ergies, it is conceivable that the parton density locally 
in the proton can become comparable to that in the nu- 
cleus. In the semi-dense/pA case, one solves the Yang- 
Mills equations [£> M , F^] — J v with the light cone sources 
ju,a = S v+ 5(x~) p%(x±) + S"- S(x+)p a A (x ± ), to deter- 
mine the gluon field produced-to lowest order in the pro- 
ton source density and to all orders in the nuclear source 
density. The inclusive gluon production cross-section, in 
this framework, was first computed by Kovchegov and 
Mueller 0J and shown to be k± factorizable in Ref. |42| . 
In Ref. 03], the gluon field produced in pA collisions was 
computed explicitly in Lorentz gauge d^A^ — 0. . Our re- 
sult is exactly equivalent to that of Dumitru & McLerran 
in A T = gauge 05] . The well known "Cronin" effect is 
obtained in our formalism and can be simply understood 
in terms of the multiple scattering of a parton from the 
projectile with those in the target. The Cronin effect and 
its evolution with rapidity will be discussed in the next 
section. 

Quark production in p/D-A collisions can be computed 
with the gauge field in Lorentz gauge 03]- The field is de- 
composed into the sum of 'regular' terms and 'singular' 
terms; the latter contain 5(x + ). The regular terms are the 
cases where a) a gluon from the proton interacts with the 
nucleus and produces a gg-pair outside, b) the gluon pro- 
duces the pair which then scatters off the nucleus. Naively, 
these would appear to be the only possibilities in the high 
energy limit where the nucleus is a Lorentz contracted 
pancake. However, in the Lorentz gauge, one has terms 
identified with the singular terms in the gauge field which 
correspond to the case where the quark pair is both pro- 
duced and re-scatters in the nucleus! 

Our result for quark pair production were discussed by 
Gelis and Fujii 05] at this conference. Related work for 
single quark production was also discussed by Tuchin |46) 5 . 
Unlike gluon production, neither quark pair-production 
nor single quark production is strictly k± factorizable. 
The pair production cross-section can however still be 
written in k± factorized form as a product of the un- 
integrated gluon distribution in the proton times a sum 
of terms with three unintegrated distributions, g , g , 4>qq,g 
and 4>qq. q q- These are respectively proportional to 2-point, 

3- point and 4-point correlators of the Wilson lines we dis- 
cussed previously. For instance, the distribution <p q q, g is 
the product of fundamental Wilson lines coupled to a qq 
pair in the amplitude and adjoint Wilson lines coupled 
to a gluon in the complex conjugate amplitude. For large 
transverse momenta or large mass pairs, the 3-point and 

4- point distributions collapse to the unintegrated gluon 
distribution, and we recover the previously discussed k±- 

5 In addition, see related work in Ref. |47|- for a recent review 
of k± factorization in heavy quark production, see Ref. |ri(i) . 



factorized result for pair production in the dilute/pp-limit. 
Single quark distributions are straightforwardly obtained 
and depend only on the 2-point quark and gluon correla- 
tors and the 3-point correlators. For Gaussian sources, as 
in the MV-model, these 2-, 3- and 4-point functions can be 
computed exactly as discussed in Ref. 02] and in Ref. 0JJ • 
The results for gluon and quark production in p/D- 
A collisions, coupled with the previous results for inclu- 
sive and diffractive '51 ,65,64.671 distributions in DIS sug- 
gest an important new paradigm. At small x in DIS and 
hadron colliders, previously interesting observables such 
as quark and gluon structure functions are no longer the 
right observables to capture the relevant physics. Instead 
they should be replaced by these dipole and multipole corre- 
lators of Wilson lines that seem ubiquitous in all high en- 
ergy processes and are similarly gauge invariant and pro- 
cess independent. The renormalization group running of 
these operators are a powerful and sensitive harbinger of 
new physics. 

4.3 Gluon and quark production in the dense/AA 
region (p A1 jk\ = pAilk\ ~ I). 

In nucleus-nucleus collisions, p\^.lk\ ~ I. There is no 
small expansion parameter and one has thus far not been 
able to compute particle production analytically in the 
CGC. Unlike gluon production in the pp and pA cases, 
fc_i_-factorization breaks down in the AA-case 69 A HHJ- k± 
factorization breaking terms are 0(1) and there are a large 
number of these. This is because the classical field comes 
in with a factor l/g-thus each insertion on the gluon is of 
order 0(1). A significant consequence is that one cannot 
factor the quantum evolution of the initial wavefunctions 
into unintegrated gluon distributions unlike the pA case. 

Nevertheless, there is a systematic way to include small 
x effects in the A A case. The problem of nuclear collisions 
is well defined in weak coupling and can be solved nu- 
merically |69ll7QII7l] . The numerical simulations thus far 
assume Gaussian initial conditions as in the MV model. 
These are good initial conditions for central Gold-Gold 
collisions at RHIC where the typical x is of order 1CP 2 . 
They are not good initial conditions at the LHC where 
the typical x at central rapidities will be at least an order 
of magnitude lower. In that case, one has to use solutions 
of JIMWLK RG equations [T2J . The numerical lattice for- 
malism of Ref. |69| is ideal for computing particle pro- 
duction in the forward light cone by matching the Wilson 
lines from each of the nuclei on the light cone. 

We restrict ourselves to discussing numerical solutions 
with Gaussian initial conditions. The saturation scale Q s 
(which is an input in the numerical solutions in this ap- 
proximation) and the nuclear radius R are the only param- 
eters in the problem. The energy and number respectively 
of gluons released in a heavy ion collision of identical nu- 
clei can therefore be simply expressed as 

1 dE _ ce_ ^3 
ttR? dr\ g 2 s ' 

ttR 2 di] g 2 ^ s ' ^ 
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where (up to 10% statistical uncertainity) we compute 
numerically ce = 0.25 and c/v = 0.3. Here 77 is the space- 
time rapidity. 

The number distributions of gluons can also be com- 
puted in this approach. Remarkably one finds that a) the 
number distribution is infrared finite, and b) the distribu- 
tion is well fit by a massive Bose-Einstein distribution for 
k±/Q s < 1-5 GeV with a "temperature" of ~ 0A7Q S and 
by the perturbative distribution Qg/kj_ for k±/Q s > 1.5. 

5 What the CGC tells us about matter 
produced in D-A and A-A collisions at RHIC 

Gluon distributions computed in pA collisions in the MV 
model exhibit the Cronin effect [501151) . One can show that 
this is exactly equivalent to the Glauber picture where 
partons from the proton acquire transverse momenta from 
multiple scattering off partons in the target |49II43| . How- 
ever, unlike the Glauber picture, quantum evolution in the 
CGC (see the discussion in Section OJl predicts that the 
Gaussian approximation breaks down completely when 
the X2 in the target is such that ln(l/a;2) ~ 1/ocs- I 11 
other words, as one produces gluons further and further 
forward in the proton fragmentation region (recall xf = 
xi — X2), the Glauber picture should break down. Indeed, 
that is precisely the trend that is observed in the RHIC 
Dcuteron-Gold experiments The rapid depletion of 
the Cronin effect is likely due to the onset of BFKL evolu- 
tion, while the subsequent saturation of this trend reflects 
the onset of saturation effects |42U27U53ll54ll43ll57j . Fur- 
ther, there is a natural explanation for the dramatic inver- 
sion of the centrality dependence that one observes in the 
RHIC data- it arises due to the onset of BFKL anoma- 
lous dimcnsions-crudely put, the nuclear Brcmsstrahlung 
spectrum changes from Q 2 s /k\ — ► Q s /k±. Finally, an ad- 
ditional piece of evidence that can be adduced in support 
of the CGC picture is the broadening of azimuthal cor- 
relations [23: for which preliminary data now exists from 
the STAR collaboration [HHj- These ideas can be tested 
conclusively in photon and di-lepton production in D-A 
collisions at RHIC H23E1E21- 

We have learnt several things from applying the CGC 
to D-A collisions. Firstly, the Gaussian MV model works 
at mid-rapidities x ~ 10 _2 -quantum evolution a la BFKL 
is not significant at these values of x. The MV model is 
therefore a good model of the initial conditions for A-A 
collisions at RHIC. (More on this to come shortly.) The 
MV model is not a good model as one goes forward in 
the Deuteron direction-at small x's of x ~ 10~ 3 or lower. 
Quantum evolution effects, seen explicitly in solutions of 
the B-K equation- are important. They will therefore pro- 
vide the initial conditions for heavy ion collisions at the 
LHC. 

The MV model when applied to heavy ion collisions 
correctly predicted the initial multiplicity at RHIC |69| . 
It was also remarkably successful in explaining rapidity 
distributions and the centrality dependence of multiplic- 
ities [23]. However, it soon became clear that the CGC 



alone was not sufficient to explain the RHIC data since a) 
it could not explain the RHIC V2 data [ZOj and b) it pre- 
dicted a suppression in D-A collisions at RHIC (the MV 
model notwithstanding) which disagreed with the RHIC 
data 56 . This failure of the CGC (here meaning quan- 
tum evolution as opposed to the MV model which has 
no evolution) thus strongly suggested that final state in- 
teractions are important at RHIC-which corroborates the 
remarkable success of hydrodynamic models. 

Why do predictions of bulk features-the multiplicity 
and rapidity and centrality dependence |73| do so well 
then? If hydrodynamic behavior sets in early, and vis- 
cous effects are small, the bulk features from the initial 
conditions will be preserved by hydrodynamic flow. This 
is seen in the hydrodynamic simulations of Hirano and 
Nara [73] • Thus one has the beginnings of a consistent 
phenomenological picture-though many puzzles remain. 
We don't understand why thermalization is early (more 
on that in the next section) or have a quantitative under- 
standing of what the viscous corrections are. In general, 
we don't have a good understanding of the properties of 
the strongly interacting quark gluon plasma- better data 
and better calculations will be helpful. 

The RHIC data on the multiplicity (approximately 
1000 hadrons in one unit of rapidity) and transverse en- 
ergy (approximately 500 GeV for central rapidities) of 
produced hadrons combined with Ea. ll2l Dlace strong con- 
straints on what Q s can be. If Q s is too small, we find, 
absurdly, that the initial transverse energy is less than the 
final measured transverse energy. If Qs is too large, we find 
that the initial multiplicity of gluons is greater than the 
final multiplicity of hadrons. While there is no obvious 
theorem that prohibits the initial gluon multiplicity being 
greater than the final hadron multiplicity, such a situation 
is unlikely in all statistical/hydrodynamic scenarios of the 
RHIC collisions. These constraints therefore allow us to 
place the bound that |72^ . 

1.3 < Q s < 2 GeV (13) 

This bound is consistent with an A 1 / 3 extrapolation of the 
Golec-Biernat-Wusthoff fit of Q 2 S to the HERA data 
A simple extrapolation gives Q s ~ 1.4 GeV. 

6 Thermalization: from CGC to QGP 

The transition to the QGP from the CGC remains as 
an outstanding theoretical problem. Due to the rapid ex- 
pansion of the system, the occupation number of modes 
falls well below one on time scales of order 1/Q S . From 
these times onwards, one expects the canonical classical 
approach to break down-well before thermalization. On 
the other hand, for elliptic flow from hydrodynamics to 
be significant, the conventional wisdom is that thermal- 
ization should set in early A necessary condition is that 
momentum distributions should be isotropic. The CGC 
initial conditions are very anisotropic with < p± >^ Q s 
and < p z >^ 0. How does this isotropization take place? 
All estimates of final state re-scattering of partons formed 
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from the melting CGC, both from 2 — ► 2 processes jZHI and 
2 — > 3 processes [JH] suggest thermalization takes longer 
than what the RHIC collisions seem to suggest-in the lat- 
ter case, Thermal ~ 13/5 g^j which at RHIC energies gives 

a s 

^thermal ~ 2 — 3 fm. 

Recently, it has been suggested that collective insta- 
bilities |77| . analogous to the well known Weibel instabil- 
ities in plasma physics, can speed up thcmalization |781 
I79| . Starting from very anisotropic (CGC-like) initial con- 
ditions, these instabilities drive the system to isotropy on 
very short time scales, of order 1/Q S in some estimates. 
What is the relation of this language of instabilities and 
that of our classical field simulations? One possibility is 
that our particular initial conditions, the non-linearities 
of the fields and the rapid expansion of the system kill 
the growth of instabilities. Another intriguing possibility 
is that small violations of boost invariance provide the the 
seeds for the instabilities. Further, to properly study ther- 
malization, one should better understand the interaction 
of high momentum (particle) and low momentum (field) 
degrees of freedom and their evolution. This leads to a real 
time renormalization group description |80| . 

An equally interesting problem is that of chemical equi- 
libration. At high energies, the initial state in a heavy ion 
collision is dominated by gluons. Are quarks produced in 
sufficient numbers for the system to reach chemical equi- 
librium (where the ratio of gluons to quarks is expected 
to be 32/21iV/)? One would expect, in weak coupling, 
that the production of quarks to be suppressed. However, 
since the fields from the CGC are of order 1/g, strong 
fields could drive the system to chemical equilibrium. First 
steps have been taken to study this problem |81U82| which 
involves numerically solving the Dirac equation in the 
background field of the two nuclei. One expects further 
progress on this problem in the near future. 



7 Open Issues in the CGC 

The CGC is a framework to think about problems in high 
energy QCD. There are many loose ends. A topic of much 
excitement among theorists recently is whether there are 
contributions beyond the JIMWLK equations-in particu- 
lar those that generate " Pomeron loops" . These contribu- 
tions are likely at low to moderate parton densities where 
impact parameter fluctuations are large. This topic is ad- 
dressed in the talk by Iancu |32]. We addressed the issue of 
fc_i_-factorization and why "dipole" and "multipole" oper- 
ators may be more relevant variables at high energies than 
structure functions. Can one derive factorization theorems 
in this framework analogous to those derived previously 
for Collinear Factorization? Turning to phenomenology, 
we have the beginnings of a consistent phenomenological 
picture of the CGC and the QGP in D-A and A-A colli- 
sions. For this to become a quantitative science, we need 
to understand the problem of thermalization from first 
principles in QCD. It is a difficult task but by no means 
an impossible one. 
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